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For the statistics of global observables in disordered systems, we discuss the matching between 
typical fluctuations and large deviations. We focus on the statistics of the ground state energy Eq 
in two types of disordered models : (i) for the directed polymer of length N in a two-dimensional 
medium, where many exact results exist (ii) for the Sherrington-Kirkpatrick spin-glass model of N 
spins, where various possibilities have been proposed. Here we stress that, besides the behavior of the 
disorder-average E§ V (N) and of the standard deviation AEq(N) ~ N"f that defines the fluctuation 
exponent Uf, it is very instructive to study the full probability distribution II(ii) of the rescaled 
variable u = E °^e^n) ^ ■ ( a ) numerically, the convergence towards IT(u) is usually very rapid, 
so that data on rather small sizes but with high statistics allow to measure the two tails exponents 
n± defined as lnn(it — > ±00) ~ —\u\ v± . In the generic case 1 < rj± < +00, this leads to explicit 
non-trivial terms in the asymptotic behaviors of the moments of the partition function when 
l/") ' the combination [|n| JV'"-'' ] becomes large (b) simple rare events arguments can usually be found to 

obtain explicit relations between n± and ujf. These rare events usually correspond to 'anomalous' 
large deviation properties of the generalized form R(w± — ^^'^ — ) ~ e - ArP±TC ±(™±) (the 
'usual' large deviations formalism corresponds to K± = 1 = p±). 

I. INTRODUCTION 

In the field of disordered systems, the interest has been first on self-averaging quantities, like the free-energy per 
degree of freedom, or other thermodynamic observables that determine the phase diagram. However, it has become 
clear over the years that a true understanding of random systems has to include the sample-to-sample fluctuations of 
global observables, in particular in disorder-dominated phases where interesting universal critical exponents show up. 
Besides these typical sample-to-sample fluctuations, it is natural to characterize also the large deviations properties, 
since rare anomalous regions are known to play a major role in various properties of random systems. 
1^1 ' Among the various global observables that are interesting, the simplest one is probably the ground-state energy Eq 
of a disordered sample. Since it is the minimal value among the energies of all possible configurations, the study of its 
CsJ ■ distribution belongs to the field of extreme value statistics. Whereas the case of independent random variables is well 
classified in three universality classes [l|, the problem for the correlated energies within a disordered sample remains 
open and has been the subject of many recent studies (see for instance 0] and references therein). For many-body 
models with N degrees of freedom (N spins for disordered spin models, N monomers for disordered polymers models), 
the interest lies 

(i) in the scaling behavior of the average Eq v (N) and the standard deviation AEq(N) with TV. Following the 
definitions of Ref. [3|, the 'shift exponent' lo s governs the correction to extensivity of the averaged value 

Eg°(N)~Neo + N u -e 1 + ... (1) 
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whereas the 'fluctuation exponent' uif governs the growth of the standard deviation 

AE (N) - N^i (2) 
(ii) in the asymptotic distribution H(u) of the rescaled variable 

_ E -E§ V (N) 



AEo(N) 



(3) 



in the limit N — > 00 



This scaling function H(u) describes the typical events where the variable u is finite. 

(iii) in the large deviations properties. In the standard 'large deviation formalism' (see for instance the recent 
review [J] and references therein), one is interested in the exponentially rare events giving rise to a finite difference v 
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between the intensive observable (Eq/N) and its averaged value Eq v (N)/N 

In disordered systems, the probability distributions of these rare events is not necessarily exponentially small in N 
but can sometimes involve other exponents 7± (see examples below in the text) 



■ X 



In this paper, we discuss these properties for two types of disordered models : for the directed polymer of length 
N in a two-dimensional medium, where many exact results exist, and for the Shcrrington-Kirkpatrick (SK) spin-glass 
model of N spins, where various possibilities have been proposed from numerical results or theoretical arguments. 
The main conclusions we draw from these two cases are the following : 

(a) it is very instructive to study the tails of the full probability distribution Tl(u) of Eq. [4] : these tails are usually 
described by the following form 

lnIT(w) oc -(-it)"- + ... 

u— > — oo 

lnIT(w) oc -u r '+ + ... (7) 

u— >+oo 

where the two tails exponents rj± are usually different and in the range 1 < r/± < +oo. In particular, the very common 
fits based on generalized Gumbel distributions are very restrictive and very misleading since they correspond to the 
unique values r\- = 1 and r/ + = +oo. We also discuss the consequences of Eq. [7] for the moments Z 1 ^ of order n 
(either positive or negative) of the partition function Zn at very low temperature. 

(b) simple rare events arguments can usually be found to obtain explicit relations between rj± and ujf. The 
probability distributions of these rare events usually correspond to 'anomalous' large deviation properties of the 
generalized forms 



_ E -Eg v {N) 

i N^ W -~ N n- <[) > Ar - oo e 



R+(w + = E ° ~ EriN) > 0) * e- NP+n +^ (8) 

NK + N K + N^oo K ' 

The paper is organized as follows. In Section m we recall the exact results concerning the directed polymer in 
a two-dimensional random medium, and discuss their meaning for the above points (a) and (b) . In Section IIII1 we 
discuss the case of the Sherrington-Kirkpatrick spin-glass model, and we present numerical results obtained for small 
sizes but with high statistics. Our conclusions are summarized in section HVl 

II. REMINDER ON THE DIRECTED POLYMER IN A TWO-DIMENSIONAL RANDOM MEDIUM 

A. Brief summary of exact results 

The directed polymer model in a two-dimensional random medium (see the review Q) is an exactly soluble model 
that has the following properties : 



(i) a single exponent 



\ (9) 



governs both the correction to extensivity of the average Efi v (N) (Eq. [TJ) and the width AE (N) (Eq. [2]). 

(ii) the rescaled distribution LT(tt) of Eq. |4] is the Tracy- Widom distribution of the largest eigenvalue of random 
matrices ensembles [8T4l0|. In particular, the two tails exponents of Eq. [7] read 

3 

V - =2 

V+ = 3 (10) 
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(iii) the exponents of the large deviations forms of Eq. [6] are respectively [lll-fli} 

7 _ = 1 

7+ =2 (11) 
After this brief reminder of known results, we now turn to their physical interpretation. 

B. Physical interpretation of the large deviation exponents in terms of simple rare events 

As explained in detail in |l4( , the large deviation exponents of Eq. [10] can be understood as follows 

(-) to obtain a ground state energy which is extensively lower than the typical, it is sufficient to draw N anomalously 

good on-site energies along the ground state path. This will happen with a probability e ~( cst ) N corresponding to 

7 _ = 1 of Eq. HI] 

(+) to obtain a ground state energy which is extensively higher than the typical, one needs to draw N 2 bad on-site 
energies (i.e. in the whole sample). This will happen with a probability e~( cst ' N corresponding to 7+ = 2 of Eq. [IT] 

Note that in the Asymmetric Exclusion process language, the interpretation is that to slow down the traffic, it is 
sufficient to slow down a single particle, whereas to speed up the traffic, one needs to speed up all particles [11| . In 
the random matrix language, the interpretation is that to push the maximal eigenvalue inside the Wigner sea, one 
needs to reorganize everything, whereas to pull the maximal eigenvalue outside the Wigner sea, one may leave the 
Wigner sea unchanged for the other eigenvalues [H, EH . 

The fact that these large deviation exponents j± can be guessed via simple physical arguments is an important 
lesson that is very useful in other disordered models which are not exactly solvable : in particular, these arguments can 
be easily extended to the directed polymer in a random medium of higher dimensionality [l4j , or to other observables 
in various models EMI. 



C. Matching between typical fluctuation and large deviations 



For an arbitrary probability distribution, the typical fluctuations in the bulk and the rare fluctuations in the far 
tails are in general different questions. However, for the probability distribution of the ground state energy Eq(N) (or 
more generally the probability distributions of other global observables) in disordered statistical physics models, it is 
very natural, from a physical point of view, to expect some matching between the typical fluctuation scaling regime 
where Eq — Eq v ~ N^f and the large deviations scaling regime where Eq — Eft v ~ TV. More precisely, the tails in the 
regions u — > ±00 of the rescaled distribution Tl(u) of typical fluctuations should match smoothly the large deviation 
region regime where the variable v of Eq. [5] is finite, which corresponds to the regime where the variable u of Eq. [3] 
is of order u oc AT 1- "' . If one plugs this scaling into the asymptotic form of Eq. [7] and if one insists that one should 
then recover the large deviations exponents of Eq. [5] one obtains the very simple relations between exponents [14j 

(l-w/h±=7± (12) 

For the directed polymer in a two-dimensional random medium, these relations are satisfied by the values quoted in 
Eqs 191101 and 1111 This smooth matching has also been discussed in the equivalent problems concer ning the current in 
the asymmetric exclusion process and the largest eigenvalue of Gaussian random matrices [T2I Il3|. 

This matching property between typical fluctuation and large deviations is again an important lesson that can be 
used in other disordered models which are not exactly solvable. These relations have been checked in detail for the 
directed polymer in dimension d — 2, 3 [l?} as well as on hierarchical lattices [l4j]. This matching property has also 
been used recently for the distribution of the dynamical barriers [IH, EH . It is also interesting from a physical point 
of view, because the asymmetry ?}_ < r/ + seen in the distribution of typical events can be seen as a consequence of 
the asymmetry of rare events 7_ < 7+. 



D. Consequences for the moments of the partition function 

Since a direct calculation of the probability distribution of the ground state of a disordered model is usually very 
difficult, analytical calculations usually focus on the moments ZZ of the partition function Z^- Then one can use 
two types of arguments to relate the distribution of Zn = e~^ FN to the distribution Pn(Eo) of the ground state 
energy Eq(N) : (1) at very low temperature T — > 0, the partition function will be dominated by the ground state 
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Zn{P +00) ~ e~P E °( N ) (2) moreover in some models, where the disorder-dominated phase < T < T c is governed 
by a zero-temperature fixed point, one expects that in the whole region of temperatures < T < T c , the probability 
distribution of the free-energy Fn will actually have the same properties as the distribution of Eq. Since (2) is valid 
for the directed polymer model, but cannot be taken for granted in all disordered models, we will restrict here to the 
point of view (1) of very low temperature T — > 0. 

There exists a simple argument that has been proposed by Zhang [5j on the specific case on the directed polymer, 
that relates the scaling behaviors of the moments with the size TV and with the replica index n to the properties 
of Pn{Eq). The idea is to evaluate the moments by using the rescaled distribution of Eq. H] 

Z% = j dE Q P N (E )e^ nE ° ~ * e-^W J duU(u)e-^ nAE ^ (13) 

For the case n > considered by Zhang Q , the integral can be then evaluated by a saddle point method in the region 
u — > — oc, where one may use the asymptotic behavior of Eq. [7] with the exponent 77 _ > 1 : the saddle point is of 
order 



u* oc {PnAEo(N))"-- 1 (14) 

that should be large u* ^> 1, and one obtains 



Z n>0 ^ e -PnE^(N) + (cst){(3n&E (N))^-- L ^ 

For the case n < 0, the equivalent calculation yields in term of the other tail exponent 77+ > 1 : 



Z n<0 ^ e ~fSnE^(N) + (cst)(f}(-n)AE (N))^+- 1 ^g-j 

For the directed polymer in a two-dimensional random medium, one obtains, using AEq(N) oc N u f with the explicit 
values of EqsHlO] 



Z™>° oc e -P nE S V (A r ) + (cst)( / 3r l )' 1 N ( 17 ) 

where one recognizes the combination (n 3 N) that appears in the Bethe Ansatz replica calculation of Ref. [7\ . Moreover, 
in Zhang's argument [5], one actually imposes that the non-trivial term of Eq. [T5]should be extensive in TV (because for 
positive integer n, the moments of the partition function can be formulated in terms of the iteration of some transfer 
matrix, and thus they have to diverge exponentially in N with some Lyapunov exponent) to obtain the relation 
u)f v Zi\ = 1 (which is equivalent here to the relation of Eq. [T^] obtained previously by the rare event interpretation). 
For n < 0, the obtained behavior 



Z™ <0 OC e -PnES v (N)+(cst)(J3(-n))-ZNl ^g) 

is rather different : the only extensive contribution of order N in the exponential comes from E§ V (N). The leading 
contribution due to fluctuations is only of subleading order TV 2, and it involves a non-integer power of the replica 
index (—n). To the best of our knowledge, the behavior of these negative moments has not been much discussed in 
the literature, in contrast to the case n > 0. 

These saddle-point calculations based on the facts that the tails exponents satisfy rj± > 1 can be very useful in 
other non-exactly soluble models, for instance in the Sherrington-Kirkpatrick spin-glass model that we now consider. 



III. SHERRINGTON-KIRKPATRICK SPIN-GLASS MODEL 



For short-ranged spin-glasses in any finite dimension d, it has been proven that the fluctuation exponent of Eq. [2] 
is exactly 9f = d/2 18]. Accordingly, the rescaled distribution IT(w) of Eq. was numerically found to be Gaussian 
in d = 2 and d = 3 [3], suggesting some Central Limit theorem. On the contrary, in mean-field spin-glasses, the width 
does not grows as TV 1 / 2 and the distribution is not Gaussian, as will be discussed in more details in this section. 
Studies on long-ranged one-dimensional spin-glasses [l9| have confirmed that non-mean-field models are characterized 
by Gaussian distributions, whereas mean-field models are not. 
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A. Brief summary of previous works 

The statistics of the ground state energy of the Sherrington- Kirkpatrick spin-glass model [2(| 

H = - J2 J v S i S o ( 19 ) 

\<i<j<N 

where the couplings J^- are random quenched variables of zero mean J = and of variance J 2 = 1/N, has been much 
studied recently with the following outputs : 

(i) there seems to be a consensus (see for instance [H, [2l|, [22| and references therein) on the shift exponent of Eq. Q] 

Us = I (20) 
whereas the 'fluctuation exponent' w/ is still under debate between the value (see [HHH [23[ and references therein) 

<"/ - \ (21) 



and the value (see [21|, 12414261 ] and references therein) 



Uf = - (22) 



(ii) the asymptotic distribution II(u) of Eq. @] has been measured numerically by various authors (see [2lj, |27|, |28[ 
and references therein), but unfortunately it has almost always been fitted by 'generalized Gumbel distributions' of 
the form e m '"~ e < containing a single free-parameter m for the shape. However these fits are very restrictive and very 
misleading since the tails exponents are fixed to be 

GeneralizedGumbel : r\- =1 

r/+ = +oo (23) 

for any value of the parameter m. In this paper, we propose instead that these exponents are in the range 1 < r]± < 
+oo. 

(iii) the large deviation properties have been also very controversial. In (29| , numerical results have been interpreted 
with the following values for the exponents y± of Eq. [5] 

Ref.[29] : 7 _ ~ 1.2 

7+ * 1.5 (24) 

Other proposals are (see [26| and references therein) 

Ref.[26] : 7_ =1 

7+ - 2 (25) 

After this brief summary of conflicting proposals, we now turn to the analysis along the same line as in the previous 
section concerning the directed polymer model. 

B. Discussion of simple rare events 

The simplest rare events one may consider for the SK model are the following : 

(-) to obtain a ground state energy which is much lower than the typical, it is natural to consider the anomalous 
ferromagnetic samples [29| that appear with a small probability of order e -( cst ) N (one needs to draw iV 2 positive 
couplings in Eq. IT5)) . and that will corresponds to anomalously low energy of order Eq (x—(cst)N 3 ^ 2 . These events 
corresponds to the 'very large deviation' of the generalized form of Eq. [5] with the values [29| 



Ferro : p_ =2 
3 



2 W 
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This form has been checked numerically in (29j . 

(+) to obtain a ground state energy which is much higher than the typical, one could consider the anomalous 

antiferromagnetic samples that appear with a small probability of order e - ^*^ 2 (one needs to draw iV 2 negative 
couplings) and that will give an energy extensively higher. In the large deviation form of Eq. [6j this would corresponds 
to 

Antiferro : 7+ = 2 (27) 

This value corresponds to the proposal of Eq. [55] from Ref. (2(| , but disagrees with the numerical proposal of Eq. [M] 
from Ref. [29]. The question is whether to obtain an extensively higher energy, it is sufficient to draw anomalously 
only a number of order N 1 - 5 random couplings instead of TV 2 . We are presently not aware of any simple argument in 
favor of this smaller power N 1 - 5 . 



C. Matching between typical fluctuation and large deviations 



In the (+) region, the matching between typical fluctuation and rare events leads to the same relation as in Eq. [12] 

(1-^)^=7+ (28) 
In particular, the possible values of 7+ and uif leads to the following values for the tail exponent ?/ + : 

If 7+ = 2 : ^+ ,= ^ = ~ = 2.6666.. 

o 

rt f=i) = — = 2.4 (29) 
5 

or 

:,<"'-» =2 

1,1"'"" = - = 1.8 (30) 
5 

In the (-) region, the matching between typical fluctuation and the very large deviations of Eq. [24] leads to the 
relation 

(k_ - Uf) 77- = p- (31) 
Using the values of Eq. [25] one obtains the two possible values for ?7_ 

if 1 8 1 a 

If u f = - -V- = 5 = 1-6 

If u f = l - r?_ =|=1.5 (32) 

If this matching works, the region of large deviation of Eq. [5] which is between the typical region and the very large 
deviation region, is constrained by consistency to involve the exponent 

7_ = (l- W/ )»?_ (33) 

The two possible values read 

If UJ = \ : 7 - =§ = 1-2 

If w/ = i : 7 - =^ = 1 -25 (34) 

Both are close to the numerical value of Eq. [M] proposed by Ref. [29( . Both disagree with the value 7_ = 1 of Eq. 
[25l used in replica calculations of [26| . 
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D. Consequences for the moments of the partition function 

As explained in detail in section III D[ the moments of the partition function Zn are then expected to follow Eqs 
ED and [16] 



Z n>0 x e -/3n_E° l, (Ar) + ( cs t)(/3„Af-/)"-- 1 
n± 

Z n<0 x e -t3nE^(N) + (cst)(f3(-n)N"f)"+- 1 (35) 

For positive n > : the two possible values of uj and of the associated tail exponent 77 _ (see Eq. 1521) correspond 
to the behaviors 
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If u f =- :Z% >0 oce-^m+Mltfrfs)^ (36) 

We note that in both cases, the non-trivial part is sub-extensive in N, in contrast to the replica calculations of 26], 
but in agreement with the replica calculations of [13, |3l|. It is also clear that the non-trivial part (n^N 1 / 2 ) for the 
case CJf = |, is simpler than the term (n 8 / 3 jV 2 / 3 ) for the case Wf = \. In both cases, the powers of n that appear are 
different from the value n 5 of perturbative replica calculations [26|, [32[ • 

For negative n < : for the case 7+ = 2 of Eq. [511 the possible behaviors are 

1 



If w/=4 -Z N 

If W/ = i -.Z^ ^-WW+MW-o)*^ (37) 

For the case 7+ =3/2 of Eq. |2"^1 proposed in Ref. [29(, the behavior of the moments can be similarly evaluated using 
Eq. 1301 Again in all cases, the non-trivial part is sub-extensive in N, as already proposed in (33j. Concerning the 
powers of (— n), the exponent (12/7) for the case u>f = | is in agreement with the replica calculations of (26|. 



E. Numerical results for small sizes and large statistics of samples 




FIG. 1: Rescaled probability distribution of the ground state energy Eo(N) in the SK model : (a) the histograms ITjv(u) of the 
rescaled variable it of Eq. [3] measured for even sizes in the range 6 < N < 20 almost coincide : this shows that the convergence 
in N towards the asymptotic form is very rapid, (b) same data in log-scale to see the tails : one sees that the left tail does not 
change, whereas finite-size effects are visible on the right tail. 



8 




FIG. 2: Tails exponents of the rescaled distribution Iljv(ii) : examples of fits of our numerical rescaled histogram (step 
function) corresponding to the size N = 10 (a) The smooth curve corresponds to the best three-parameter fit in the range 
u < — 1 by the form (ao — ai(— u) v ~) : the left tail exponent is of order r/- ~ 1.5. (b) The smooth curve corresponds to the 
best three-parameter fit in the range u > 1 by the form (ao — aiu v+ ) : the right tail exponent is of order 77+ ~ 2.8. 



Most numerical works on the distribution of the ground state energy in the SK model have followed the strategy to 
study the biggest sizes N as possible, to measure the ave rage d value and the variance (see 0, HH HH and references 
therein). An opposite strategy has been followed in Ref. [27| where an exact enumeration of the disordered samples 
with the binomial distribution ±J was performed for small sizes. As mentioned in f22j , the results for the rescaled 
distribution H(u) at N — 9 are already very good when compared to the results for larger N. In other cases, we have 
also found that the distribution of rescaled variables converge much more rapidly than other observables fla . HH l34j . 
In the following, we thus follow the same strategy : we study the distribution of Eq(N) for the small sizes with a high 
statistics of disordered samples. 

On Fig. [1] (a), we show the measured histograms IIjv(u) of the the rescaled variable u of Eq. [3] for even sizes in 
the range 6 < N < 20 with a statistics of 2.10 8 > n s (N) > 13. 10 4 disordered samples : one clearly sees that all these 
histograms almost coincide. Our conclusion is thus that the convergence in N towards the asymptotic form is very 
rapid, so that these small-size data should provide a reliable measure of the asymptotic II(u). As explained before, 
we are mainly interested into the tails exponents rj± of Eq. [7] : as shown on Fig. [T] (b) the convergence of the left tail 
is extremely good, whereas the convergence of the right tail presents much stronger finite-size effects. 

Let us first consider the left tail. The three-parameter fit of II(u) in the range u < — 1 by the form (ao — a\{— u) v -) 
yields the value (see Fig [5] (a)) 

r]- ~ 1.5 (38) 

that corresponds exactly to the value associated to w/ = 1/6 (see Eq. [32]) . Of course, it is probably not far enough 
from the alternative value ?7_ =1.6 corresponding to w/ = 1/4 (see Eq. I3"2l to really rule out the value uif — 1/4. 

Let us now turn to the right tail. The three-parameter fit of II(m) in the range u > 1 by the form (ao — aiu n+ ) yields 
values for 77+ that are less precise, as a consequence of the finite-size effects visible on Fig. Q] (b). We have already 
found in other studies that the right tail is usually more difficult to measure than the left tail [1J|. Nevertheless our 
non-precise values of rj + in the range [2.4, 2.9] (see Fig[5] (b)) seem more compatible with the value 7+ = 2 than with 
the value 7+ = 1.5 (see Eq. [2^1and [3Tj|l . 



F. Final discussion 



In summar y, e ven if a definitive agreement on the precise value of the fluctuation exponent ujf remains difficult to 
reach (see [3l l2ll |22| and references therein), our conclusions concerning the SK model are the following : 

(i) the numerical measure of the left tail exponent ?y_ is in agreement with the matching argument based on rare 
ferromagnetic samples described by the very-large deviation form of Eq. [S] with the values of Eq. [5E] from Ref. |29j | . 
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Then the large deviation form of Eq. [6] is constrained to involve an exponent 7_ given by Eq. [33] 

7 =(1-0,,),- = ^ 
(2 ~ u f> 

This explains the numerical value of Eq. [M] proposed in Ref . [2!| , and excludes the value 7_ = 1 of Eq. [55] used in 
the replica calculations of (26[. We note moreover that this 'usual large deviation value' 7_ = 1 would be satisfied 
only for the value ujf = 1/3, i.e. only if the fluctuation exponent ujf would coincide with the shift exponent u a = 1/3 
of Eq. HO] 

(ii) although less precise, the numerical measure of the right tail exponent r/ + is more in favor of the large deviation 
exponent 7+ = 2, that can be justified with a simple rare events argument (see Eq. 1271) . 

(iii) finally, the facts that the tails exponents satisfy rj± > 1 induces non-trivial behavior for the moments of the 
partition function (see Eqs. I3l)|) when (|n|iV^) becomes large. In particular, from Eqsl3"u1andl3"TI our conclusion is 
that the only extensive term in N comes from the trivial term E§ V (N) both for negative and positive n. Moreover, 
the non-trivial sub-extensive terms can a priori involve non-integer powers of the replica index n. 



IV. CONCLUSION 



In this paper, we have discussed the statistics of the ground state energy Eq(N) in two types of disordered models: 
(i) for the directed polymer of length N in a two-dimensional medium, where many exact results exist (ii) for the 
Sherrington-Kirkpatrick spin-glass model, where various possibilities are still under debate both numerically and 
theoretically. Our main conclusions are the following. Besides the behavior of the disorder-average Eq v (N) and of 
the standard deviation AE (N) ~ N u f , it is very instructive to study the full probability distribution H(u) of the 
rescaled variable u = E °^e^n) N ^ ' 

(a) numerically, the convergence towards H(u) is usually very rapid, so that data on rather small sizes but with high 
statistics allow to measure the tails exponents r)± that satisfy generically 1 < r/± < +00 (whereas the very common 
fits based on generalized Gumbel distributions correspond to the unique values 77 = 1 and ry + = +00). Moreover, if 
one wishes to measure tails beyond the region probed via simple sampling, one may uses a Monte-Carlo procedure in 
the disorder, as done in 28] for the SK model, and in [17J for the directed polymer model. 

(b) simple rare events arguments can usually be found to obtain explicit relations between r)± and w/. These rare 
events usually correspond to 'anomalous' large deviation properties of the generalized form Rn(w± — Eo ^ N ^t — ) ~ 
e -jv p ±TC ± (to ± ) > uSUa p> l ar ge deviations formalism corresponding to K± = 1 = p± is too restrictive for disordered 
models, as shown on explicit examples in the text). 

(c) We have also discussed the consequences of 1 < rj± < +00 for the moments Z% of order n (either positive or 
negative) of the partition function Zjq. In the regime where [|tt.| A'"^^ ] becomes large, a saddle-point calculation leads 
to explicit non-trivial terms in the asymptotic behaviors of the moments Z^ of the partition function. 

We have shown in detail how this analysis for the directed polymer is in agreement with all known exact results. 
For the SK model, we have explained how this analysis agrees or disagrees with various possibilities debated in the 
literature. 

In conclusion, we believe that this type of analysis based on the matching between typical fluctuations and rare 
events is very useful to study disordered systems. Here we have focused on the statistics of the ground state energy, 
but it can also be used for other global observables such as the maximal dynamical barrier of a disordered sample 
[HI EB] ! or for the statistics of large excitations in ferromagnets and spin-glasses [l4[ . 
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